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We present results of a study of neutrino oscillation based on a 766 ton-year exposure of KamLAND to 
reactor anti-neutrinos. We observe 258 Fe candidate events with energies above 3.4 MeV compared to 365.2 
events expected in the absence of neutrino oscillation. Accounting for 17.8 expected background events, the 
statistical significance for reactor Fe disappearance is 99.998%. The observed energy spectrum disagrees with 
the expected spectral shape in the absence of neutrino oscillation at 99.6% significance and prefers the dis- 
tortion expected from 77^ oscillation effects. A two-neutrino oscillation analysis of the KamLAND data gives 
^ eV'^. A global analysis of data from KamLAND and solar neutrino experiments yields 



Am^= 7.9l°:^xl0 
Am2 = 7.9t;|!jxl0- 



' eV^ and tan'^ ^=0.40^007, the most precise determination to date. 



PACS numbers: 14.60.Pq, 26.65.+t, 28.50.Hw 



The first measurement of reactor anti-neutrino disappear- 
ance by KamLAND 1 1 ] suggested that solar neutrino flavor 
transformation through the Mikheyev-Smirnov-Wolfenstein 
(MSW) 1 2] matter effect has a direct correspondence to anti- 
neutrino oscillation in vacuum. Assuming CPT invariance, 
KamLAND and solar-neutrino experiments have restricted 
the solar oscillation parameters, eliminating all but the large- 
mixing-angle (LMA-MSW) solution. This Letter reports 
more stringent constraints on neutrino oscillation parameters 
from KamLAND based on a three times longer exposure and 
a 33% increase in the usable fiducial volume. Large variations 
in the reactor power production in Japan in 2003 allowed us 
to study the V^. flux dependence. The first evidence for spec- 
tral distortion in the V^ spectrum is provided here; spectral 



distortion is direct evidence of an oscillation effect. 

KamLAND consists of 1 kton of ultra-pure liquid scintil- 
lator (LS) contained in a 13-m-diameter transparent nylon- 
based balloon suspended in non-scintillating oil. The balloon 
is surrounded by 1 879 photomultiplier tubes (PMTs) mounted 
on the inner surface of an 18-m-diameter spherical stainless- 
steel vessel. Electron anti-neutrinos are detected via inverse 
/3-decay, v^+p ^ e^ +n, with a 1.8 MeV Ve energy thresh- 
old. The prompt scintillation light from the e+ gives an es- 
timate of the incident v^ energy, Ej:^ = i^prompt + E„ + 
0.8 MeV, where E'piompt is the prompt event energy includ- 
ing the positron kinetic energy and the annihilation energy, 
and En is the average neutron recoil energy. The ~ 200 us 
delayed 2.2 MeV 7 ray from neutron capture on hydrogen is 



a powerful tool for reducing background. On average, neu- 
trons are captured within 9 cm and the spatial correlation be- 
tween prompt and delayed signals is dominated by the ver- 
tex position resolution. A 3.2kton water-Cherenkov detec- 
tor surrounds the containment sphere, absorbing 7 rays and 
neutrons from the surrounding rock and tagging cosmic -ray 
muons. This outer detector (OD) is over 92% efficient for 
muons passing through the fiducial volume. 

KamLAND is surrounded by 53 Japanese power reac- 
tor units. The reactor operation data, including thermal 
power generation, fuel burn up, exchange and enrichment 
records, are provided by all Japanese power reactors and 
are used to calculate fission rates of each isotope. The 
averaged relative fission yields for the run period were 

235u . 238u : 239p^ . 241p^ ^ q 553 . q Q79 . q jqj . q Q57 j^^^ 

expected Fg flux is calculated using fission rates and Ve spec- 
tra; the spectra are from Ref. t3J. The Ve contribution from 
Japanese research reactors and all reactors outside of Japan is 
4.5%. We assume that these reactors have the same average 
fuel composition as the Japanese power reactors. The total 
integrated thermal power flux over the detector livetime was 
701 Joule/cm^. 

We report on data collected between March 9, 2002 and 
January 11, 2004, including re-analysis of the data used in 
Ref. 1 1 ] . The central detector PMT array was upgraded on 
February 27, 2003 by commissioning 554 20-inch PMTs, in- 
creasing the photo-cathode coverage from 22% to 34% and 
improving the energy resolution from 7.3%/ y/E{MeY) to 
6.2%/v/^(MeV). The trigger threshold of 200 hit 17-inch 
PMTs corresponds to about 0.7 MeV at the detector center and 
has an efficiency close to 100%. We use a prompt event en- 
ergy analysis threshold of 2.6 MeV to avoid backgrounds in- 
cluding the effect of anti-neutrinos from uranium and thorium 
decaying in the Earth (geo-neutrinos). 

The location of interactions inside the detector is deter- 
mined from PMT hit timing; the energy is obtained from the 
number of observed photo-electrons after correcting for posi- 
tion and gain variations. Position and time dependence of the 
energy estimation are monitored periodically with 7-ray and 
neutron sources along the central vertical axis (z-axis) of the 
scintillator volume. Trace radio-isotopes on the balloon and 
in the scintillator are also exploited. The systematic uncer- 
tainty in the energy scale at the 2.6 MeV prompt event energy 
(E-jj^ ~ 3.4 MeV) analysis threshold is 2.0%, corresponding 
to a 2.3% uncertainty in the number of events in an unoscil- 
lated reactor Ve spectrum. 

The radial fiducial volume cut is relaxed from 5 m 1 1 ] to 
5.5 m in the present analysis, expanding the fiducial mass to 
543.7 tons (4.61 x 10"^^ free target protons). The radial posi- 
tions of the prompt and delayed event are both required to 
be less than 5.5 m. The 1.2 m cylindrical cut along the z- 
axis previously used to exclude low energy backgrounds from 
thermometers is not applied. The event selection cuts for the 
time difference (AT) and position difference (Ai?) between 
the positron and delayed neutron are 0.5 /is < AT < 1000/is 
and AR < 2 m, respectively. The event energies are re- 



quired to be 1.8 MeV < ii^deiayed < 2.6 MeV and 2.6 MeV < 
^prompt < 8.5 MeV. The efficiency of all cuts is (89.8 ± 1 .5)%. 

The total volume of the liquid scintillator is 1171 ± 25 m^, 
as measured by flow meters during detector filling. The 
nominal 5.5-m-radius fiducial volume (^ttR^) corresponds to 
0.595 ±0.013 of the total LS volume. The effective fiducial 
volume is defined by the cuts on the radial positions of the 
reconstructed event vertices. At present, only z-axis calibra- 
tions are available, so we assess the systematic uncertainty 
in the total fiducial volume by studying uniformly-distributed 
muon spallation products, identified as delayed coincidences 
following muons. We measure the position distribution of 
the /3-decays of ^^B (Q= 13.4MeV, T1/2 = 20.2ms) and ^^N 
((5= 17.3 MeV, ri/2= 11-Oms), which are produced at the 
rate of about 60 ^^B/^^N events/kton-day. Fits to the en- 
ergy distribution of these events indicate that the sample is 
mostly ^^B; the relative contribution of ^^N is only ^1%. 
The number of ^^B/^^N events reconstructed in the fiducial 
volume compared to the total number in the entire LS vol- 
ume is 0.607 ± 0.006(stat) ± 0.006(syst), where the system- 
atic error arises from events near the balloon surface that de- 
posit a fraction of their energy outside the LS. As a consis- 
tency check, in a similar study of spallation neutrons, which 
we identify via the 2.2 MeV capture 7 ray, we find the ratio 
0.587 ±0.013(stat). 

The ^^B/^^N events typically have higher energy than 
Ve candidates, so an additional systematic error accounts for 
possible dependence of effective fiducial volume on energy. 
We constrain the variation to 2.7% by comparing the prompt 
and delayed event positions of delayed-neutron /3-decays of 
^Li (Q= 13.6 MeV, Ti/2 = 178 ms) and ^He (Q= 10.7 MeV, 
Ti/2= 119ms). The observed capture distance variation is a 
measure of the energy uniformity of the vertex finding algo- 
rithm. Combining the errors from the LS volume measure- 
ments, the ^^B/^^N volume ratio calibration, and the con- 
straints on energy dependence, we obtain a 4.7% systematic 
error on the fiducial volume. 

The rate of accidental coincidences increases in the outer 
region of the fiducial volume, since most background sources 
are external to the liquid scintillator This background is esti- 
mated with a 10 ms to 20 s delayed-coincidence window and 
by pairing random singles events. These consistent methods 
predict 2.69 ± 0.02 events above the 2.6 MeV threshold. 

Above 2.6 MeV, neutrons and long-lived delayed-neutron 
/3-emitters are sources of correlated backgrounds. The ^3000 
spallation-produced neutrons per kton-day are effectively 
eliminated with a 2 ms veto of the entire detector following 
a detected muon. The remaining fast neutrons come from 
muons missed by the OD or interacting in the rock just outside 
it. This background is reduced significantly by the OD and 
several layers of absorbers: the OD itself, the 2.5 m of non- 
scintillating oil surrounding the LS, and the 1 m of LS outside 
the fiducial volume. We estimate this background contributes 
fewer than 0.89 events to the data sample. 

The uncorrelated background from ^^B/^^N spallation 
products is effectively suppressed by the delayed-coincidence 



TABLE I: Estimated systematic uncertainties (%). 



iducial Volume 


4.7 


Reactor power 2. 1 


nergy threshold 


2.3 


Fuel composition 1.0 


fficiency of cuts 


1.6 


Vc spectra [3] 2.5 


ivetime 


0.06 


Cross section [5] 0.2 



Total systematic uncertainty 



6.5 



requirement. However, the ^1.5 events/kton-day in the 
delayed-neutron branches of ^Li and ^He mimic the Fe sig- 
nal. From fits to the decay-time and /3-energy spectra we 
see mostly ^Li decays; the contribution of ^He relative to 
^Li is less than 15% at 90% C.L. For isolated, well-tracked 
muons passing through the detector, we apply a 2 s veto within 
a 3 m radius cylinder around the track. We veto the entire 
volume for 2 s after one in ^30 muons, those that produce 
more than ^10^ photo-electrons above minimum ionization 
or muons tracked with poor reliability. We estimate that 
4.8 ± 0.9 ^Li/®He events remain after the cuts. The dead- 
time introduced by all muon cuts is 9.7%; the total livetime 
including spallation cuts is 515.1 days. 

A third source of correlated background comes indirectly 
from the a decays of the radon daughter ^^°Po in the liquid 
scintillator The signal of the 5.3 MeV a-particle is quenched 
below the threshold, but the secondary reaction ^•^C(Q;,n)^^0 
produces events above 2.6 MeV. Special runs to observe the 
decay of ^^OPo establish that there were (1.47±0.20)x 10^ a 
decays during the livetime of data taking. Using the ^^C(a,n) 
reaction cross sections from Ref. |4], Monte Carlo simulations 
and detailed studies of quenching effects to convert the outgo- 
ing neutron energy spectrum into a visible energy spectrum, 
we expect 10.3 ±7.1 events above 2.6 MeV. The spectrum is 
concentrated around 6 MeV, from decays of levels in ^^O, and 
4.4 MeV, from 7 decays following neutron inelastic scatter- 
ing on ^^C. The observed energy from neutron-proton elastic 
scattering is mostly quenched below 2.6 MeV. This a-induced 
background was not considered in Ref. |jj] and would have 
contributed 1.9 ± 1.3 additional background events (2.8 ±1.7 
total background events). The total background to the V^- 
signal above 2.6 MeV in the present analysis is 17.8 ±7.3 
events, where the bound on the fast neutron background is 
accounted for in the uncertainty. 

In the absence of anti-neutrino disappearance, we expect to 
observe 365.2 ± 23. 7(syst) Fg events above 2.6 MeV, where 
the systematic uncertainty is detailed in Table U We ob- 
serve 258 events, confirming 77e disappearance at the 99.998% 
significance level. The average Ve survival probability is 
0.658 ± 0.044(stat) ± 0.047(syst), where the background er- 
ror has been included in the systematic uncertainty. The 
effective baseline varies with power output of the reactor 
sources involved, so the survival probabilities for different pe- 
riods are not directly comparable. Applying the new analy- 
sis on the previously reported data (March 2002 to October 
2002) IIJ gives 0.601 ±0.069(stat)±0.042(syst), in agree- 
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FIG. 1: (a) Estimated time variation of the reactor V^ flux at 
KamLAND assuming no anti-neutrino oscillation, (b) Observed 77^ 
event rate versus no-oscillation reactor V^ flux. Data points corre- 
spond to intervals of approximately equal 77^ flux. The dashed line is 
a fit, the 90% C.L. is shown in gray. The solid line is a fit constrained 
to the expected background. The reactor distance distribution for Fe 
events in the absence of oscillations is shown in the inset. 



ment with 0.589 ± 0.085(stat) ± 0.042(syst), after correction 
for the (a,n) background. 

After September 2002, a number of Japanese nuclear re- 
actors were off, as indicated in Fig. QJ. This decreased the 
expected no-oscillation 77^ flux by more than a factor of two. 
In Fig. ^ the signal counts are plotted in bins of approxi- 
mately equal v^ flux corresponding to total reactor power. For 
Am^ and tan^ B determined below and the known distribu- 
tions of reactor power level and distance, the expected oscil- 
lated Fe rate is well approximated by a straight line. The slope 
can be interpreted as the reactor-coiTelated signal and the in- 
tercept as the reactor-independent constant background rate. 
Fig.fQ) shows the linear fit and its 90% C.L. region. The inter- 
cept is consistent with known backgrounds, but substantially 
larger backgrounds cannot be excluded; hence this fit does 
not usefully constrain speculative sources of anti-neutrinos 
such as a nuclear reactor at the Earth's core |6]. The pre- 
dicted KamLAND rate for typical 3 TW geo-reactor scenarios 
is comparable to the expected 17.8 ±7.3 event background 
and would have minimal impact on the analysis of the reac- 
tor power dependence signal. In the following we consider 
contributions only from known anti-neutrino sources. 

Fig. 12^ shows the correlation of the prompt and delayed 




FIG. 2: (a) The correlation between the prompt and delayed event 
energies after cuts. The three events with -Edeiaycd ~ 5MeV are 
consistent with neutron capture on carbon, (b) Prompt event energy 
spectrum of F^ candidate events with associated background spectra. 
The shaded band indicates the systematic error in the best-fit reactor 
spectrum above 2.6 Me V. 



event energy after all selection cuts except for the i^deiayed 
cut. The prompt energy spectrum above 2.6 MeV is shown in 
Fig-Eb- The data evaluation method with an unbinned max- 
imum likelihood fit to two-flavor neutrino oscillation is sim- 
ilar to the method used previously 1 1 1 . In the present analy- 
sis, we account for the ^Li, accidental and the ^■^C(Q;,n)^^0 
background rates. For the (a,n) background, the contri- 
bution around 6 MeV is allowed to float because of uncer- 
tainty in the cross section, while the contributions around 
2.6 MeV and 4.4 MeV are constrained to within 32% of the 
estimated rate. We allow for a 10% energy scale uncer- 
tainty for the 2.6 MeV contribution due to neutron quench- 
ing uncertainty. The best-fit spectrum together with the back- 
grounds is shown in Fig.|25; the best-fit for the rate-and-shape 
analysis is Am^ = T.QtgJ x 10"^ eV^ and tan^ e = 0.46, with 
a large uncertainty on tan^ 9. A shape-only analysis gives 



Ar 



= (8.0 ± 0.5)x 10-5 eV^ and tan^ 9 = 0.76. 



Taking account of the backgrounds, the Baker-Cousins x^ 
for the best-fit is 13.1 (11 DOF). To test the goodness-of-fit 
we follow the statistical techniques in Ref. ^TJ. First, the 
data are fit to a hypothesis to find the best-fit parameters. 
Next, we bin the energy spectrum of the data into 20 equal- 
probability bins and calculate the Pearson x^ statistic (xi) 
for the data. Based on the particular hypothesis 10,000 spec- 
tra were generated using the parameters obtained from the 
data and Xp was determined for each spectrum. The con- 
fidence level of the data is the fraction of simulated spectra 
with a higher Xp- For the best-fit oscillation parameters and 
the a priori choice of 20 bins, the goodness-of-fit is 11.1% 




2.6 MeV prompt , KamLAND data 

analysis threshold best-fit oscillation 

i best-fit decay 
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L„/E_ (kmAleV) 

FIG. 3: Ratio of the observed Ve spectrum to the expectation for no- 
oscillation versus Lq/E. The curves show the expectation for the best- 
fit oscillation, best-fit decay and best-fit decoherence models taking 
into account the individual time-dependent flux variations of all re- 
actors and detector effects. The data points and models are plotted 
with Lo=180 km, as if all anti-neutrinos detected in KamLAND were 
due to a single reactor at this distance. 



with Xp/DOF = 24.2/17. The goodness-of-fit of the scaled no- 
oscillation spectrum where the normalization was fit to the 
data is 0.4% (Xp/DOF = 37.3/18). We note that the Xp and 
goodness-of-fit results are sensitive to the choice of binning. 

To illustrate oscillatory behavior of the data, we plot in 
Fig. |3] the Lg/E distribution, where the data and the best- 
fit spectra are divided by the expected no-oscillation spec- 
trum. Two alternative hypotheses for neutrino disappear- 
ance, neutrino decay |8] and decoherence |9], give dif- 
ferent Lq/E dependences. As in the oscillation analy- 
sis, we survey the parameter spaces and find the best-fit 
points at (sin^ 0, m/cT) = (1.0,0.011 Me V/km) for decay and 
(8111^261, 7°) = (1.0, 0.030 Me V/km) for decoherence, using 
the notation of the references. Applying the goodness-of-fit 
procedure described above, we find that decay has a goodness- 
of-fit of only 0.7% (Xp/DOF= 35.8/17), while decoherence 
has a goodness-of-fit of 1.8% (Xp/DOF= 32.2/17). We note 
that, while the present best-fit neutrino decay point has already 
been ruled out by solar neutrino data lioll and observation of 
SN1987A, the decay model is used here as an example of a 
scenario resulting in a 17^ deficit. If we do not assume CPT 
invariance and allow the range 0.5 < sin^ 9 < 0.75, then the 
decay scenario considered here can avoid conflict with solar 
neutrino LMI and SN1987A data 111.]. 

The allowed region contours in Am^-tan^ 9 parameter 
space derived from the Ax^ values (e.g., Ax^ < 5.99 for 95% 
C.L.) are shown in Fig.|^. The best-fit point is in the region 
commonly characterized as LMA I. Maximal mixing for val- 
ues of Am^ consistent with LMA I is allowed at the 62.1% 
C.L. Due to distortions in the spectrum, the LMA II region 
(at ATO2-2xlO-4eV2) is disfavored at the 98.0% C.L., as 
are larger values of Am^ previously allowed by KamLAND. 
The allowed region at lower Arn^ is disfavored at the 97.5% 
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FIG. 4: (a) Neutrino oscillation parameter allowed region from KamLAND anti-neutrino data (shaded regions) and solar neutrino experiments 
(lines) [12]. (b) Result of a combined two-neutrino oscillation analysis of KamLAND and the observed solar neutrino fluxes under the 



assumption of CPT invariance. 
range. 



The fit gives Am =7.9 



+0.6 



xlO" 



■ Q = 0.40^ 



including the allowed 1-sigma parameter 



C.L., but this region is not consistent with the LMA region 
determined from solar neutrino experiments assuming CPT 
invariance. 

A two-flavor analysis of the KamLAND data and the ob- 
served solar neutrino fluxes [13], with the assumption of CPT 
invariance, restricts the allowed Am^-tan^ d parameters in 
Fig. 4b. The sensitivity in Arn^ is dominated by the observed 
distortion in the KamLAND spectrum, while solar neutrino 
data provide the best constraint on d. The combined analysis 
gives Ato^ = 7.9 Ij; j X 10"^ eV^ and tan^ d = OAOtol^ 

The conclusion that the LMA II region is excluded is 
strengthened by the present result. The observed distortion 
of the spectral shape supports the conclusion that the obser- 
vation of reactor 77e disappearance is due to neutrino oscilla- 
tion. Statistical uncertainties in the KamLAND data are now 
on the same level as systematic uncertainties. Current efforts 
to perform full-volume source calibrations and a reevaluation 
of reactor power uncertainties should reduce the systematic 
uncertainties. 
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